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CONSTRAINING GRB INITIAL LORENTZ FACTOR WITH THE AFTERGLOW ONSET FEATURE AND 

DISCOVERY OF A TIGHT T - E 7 j SO CORRELATION 
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ABSTRACT 

The onset of GRB afterglow is characterized by a smooth bump in the early afterglow lightcurve 
caused by the deceleration of the gamma-ray burst (GRB) fireball by the circumburst medium. 
We make an extensive search for such a deceleration feature, either from the literature for optical 
lightcurves, or from the X-ray afterglow lightcurve catalog established with the Swift/XRT. Twenty 
optically selected GRBs and 12 X-ray selected GRBs are found to show the onset signature, among 
which 17 optically selected GRBs and 2 X-ray-selected GRBs have redshift measurements. We study 
the optical z-known sample by fitting the lightcurves with a smooth broken power-law and measure the 
width (w), rising timescale (t r ), and decaying timescale (td) at full-width-at-half-maximum (FWHM). 
Strong mutual correlations among these timescales and with the peak time (t p ) are found. The optical 
peak luminosity (L Pi o) at the lightcurve bump is anti-correlated with t p and correlated with w, indi- 
cating a dimmer and broader bump at a later peak time. The ratio t r /td is almost universal among 
bursts, but the ratio t r /t p varies from 0.3 ~ 1. The isotropic gamma-ray energy (i? 7 .i SO ) is tightly 
correlated with £ p o and t p in the burst frame. Assuming that the bumps signal the deceleration of 
the GRB fireballs in a constant density medium, we calculate the initial Lorentz factor (ro) and the 
deceleration radius (Rd) of the GRBs in the optical-selected sample. It is found that Tq are typically 
a few hundreds, and the typical deceleration radius is i?dec ~ 10 17 cm. More intriguingly, a tight 
correlation between the initial Lorentz factor and the isotropic gamma-ray energy is found, namely 
Fo ~ 195E® ?J 52 (satisfied for both the optical and X-ray z-known samples). This correlation is help- 
ful to understand GRB physics, and may serve as an indicator of Tq for other long GRBs. We find 
that the early bright X-rays are usually dominated by a different component from the external shock 
emission, but occasionally (for one case) an achromatic deceleration feature is observed. Components 
in X-rays would contribute to the diversity of the observed X-ray lightcurves. 
Subject headings: radiation mechanisms: non-thermal: gamma-rays: bursts 



1. INTRODUCTION 

The fireball model has been extensively employed 
to explain the gamma-ray burst (GRB) phenomenon 
(Meszaros 2002; Zhang & Meszaros 2004; Piran 2004). 
The observed prompt gamma-ray emission is explained 
by synchrotron (or inverse Compton) emission from the 
internal shocks in an erratic, unsteady, relativistic fireball 
(Rees & Meszaros 1994), and broadband afterglow emis- 
sion is attributed to synchrotron emission from the ex- 
ternal shock when the fireball is decelerated by a circum- 
burst medium (Meszaros & Rees 1997; Sari et al. 1998). 
In order to avoid the "compactness problem" of high 
energy non-thermal photons detected from GRBs, the 
GRB fireball is required to move relativistically towards 
earth. After the initial radiation-dominated acceleration 
phase, the fireball enters the matter-dominated "coast- 
ing" phase. The fireball keeps an approximately same 
Lorentz factor until it sweeps up a considerable amount 
mass from the ambient medium at the so-called decel- 
eration radius, after which T decreases with R (and the 
observer time t) as a power law. The factor Lorentz fac- 
tor during the coasting phase is called the initial Lorentz 
factor (Fo), which is a crucial parameter to understand 
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GRB physics, but is poorly known for most GRBs. 

Three methods have been proposed to measure the ini- 
tial Lorentz factor To. The first method is to apply the 
"compactness" argument to use the high energy photon 
cutoff energy in the GRB prompt emission spectrum to 
estimate Tq. Since so far no clear cutoff feature is ob- 
served, the common practice is to use the observed maxi- 
mum photon energy to set a lower limit on Tq (Fenimorc 
et al. 1993; Woods & Loeb 1995; Baring & Harding 
1997; Lithwick & Sari 2001). This method suffers sev- 
eral uncertainties. First, the cutoff energy depends on 
both Tq and the emission radius i? 7 (Gupta & Zhang 
2008). The method therefore relies on the assumption of 
i? 7 = T^cSt, the internal shock radius. Such an assump- 
tion is not necessarily correct. Second, the minimum 
variability time scale St is subject to large uncertainty, 
since the GRB lightcurves are chaotic and do not have 
a characteristic time scale. Finally, Fermi observations 
(Abdo et al. 2009b, c) indicate that some GRBs have a 
distinct high energy component in the GeV range, which 
may come from a different emission region. A straight- 
forward usage of the method may lead to erroneous con- 
clusions. The second method to measure Tq is to use 
the blackbody component detected in some GRB spec- 
tra (Pe'er et al. 2007). The limitation of the method is 
the difficulty of identifying blackbody components from 
the GRB spectra. So far the only case that this method 
is securely applied is GRB 090902B (Ryde et al. 2009). 
The third, but most commonly adopted, method is to use 
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the early afterglow lightcurve. The fireball model pre- 
dicts that at the onset of deceleration, the emission from 
the forward shock would display a smooth bump in the 
lightcurve, the peak of which corresponds to the epoch 
when essentially half of the fireball energy is transferred 
to the medium (Sari & Piran 1999; Kobayashi & Zhang 
2007). The most relevant case is the thin shell regime, 
defined by that the thickness of the fireball shell satisfies 
A < (^/nm p c 2 ) 1//3 r o 8 ^ 3 , where E is the kinetic energy 
of the fireball, n is the medium density, m p the mass of 
proton, c speed of light (Kobayashi et al. 1999; Sari & Pi- 
ran 1999). Within this regime, the deceleration time (the 

peak time at the lightcurve bump, t p oc T 8 ^ 3 (_E/n) 1 / 3 , 
Meszaros & Rees 1993) is sensitive depends on the initial 
Lorentz factor but insensitive to other parameters. The 
detection of this time can be then used to infer To. In the 
optical band, early emission may be contaminated by the 
emission from the reverse shock (Meszaros & Rees 1997; 
Sari & Piran 1999; Kobayashi 2000; Zhang et al. 2003). 
However, under certain conditions (either a Poynting flux 
dominated flow, Zhang & Kobayashi 2005; or a relatively 
low typical synchrotron frequency in the reverse shock, 
Jin & Fan 2007), the reverse shock component would not 
show up in the optical band. In these bursts, a smooth 
onset bump can be detected, which signals the deceler- 
ation feature of the fireball, and hence, can be used to 
constrain the initial Lorentz factor and the deceleration 
radius (Sari Piran 1999, Zhang et al. 2003; Molinari et 
al. 2007; Xue et al. 2009; Zou & Piran 2009). 

In this paper, we constrain Tq with the early GRB af- 
terglows that show the deceleration signature, and inves- 
tigate the possible correlations among deceleration pa- 
rameters (including Tq) as well as the prompt gamma- 
ray emission properties. We make an extensive search 
for the onset of afterglow signature in the optical and 
X-ray lightcurves. Our sample selection is presented in 
Section 2. The temporal characteristics and their correla- 
tions are presented in Section 3, and the relation between 
the prompt gamma-ray properties and the deceleration 
properties are investigated in Section 4. In particular, 
we constrain Tq and the deceleration radius of the fire- 
ball for the z-known sample, and discover a tight corre- 
lation between r and -E 7 ,i S o- Discussion and conclusions 
are presented in Sections 5 and 6, respectively. A con- 
cordance cosmology with parameters Hq = 71 km s _1 
Mpc -1 , Qm = 0.30, and £1a = 0.70 are adopted. Nota- 
tion Q n denotes Q/10 n in the cgs units throughout the 
paper. 

2. SAMPLE SELECTION AND LIGHTCURVE FITTING 

We make an extensive search for the smooth "bump" 
feature as the onset of GRB afterglow. The flares or the 
reversed shock emission component introduce confusion 
to identify the onset feature, especially in the X-ray band. 
We therefore employ the following two criteria to selec- 
tion the samples. First, the bump is smooth without su- 
perimposing significant flares around the bump. Second, 
the decay slopes post the peak time is shallower than —2. 
Through literature search, we obtain 20 optical-selected 
GRBs that show the afterglow onset feature. We also go 
through the Swift XRT lightcurve archive that has been 
processed by our group in the past (from January 2005 to 
September 2009), and identify 12 X-ray selected sample 



with afterglow onset signature. The observational results 
of these bursts are summarized in Tables 1 and 2. The 
details of XRT data reduction have been presented in 
a series of papers published by our group (Zhang et al. 
2007a, Paper I; Liang et al. 2007, 2008, 2009; Papers II, 
III, IV). The prompt gamma-ray properties of the bursts 
arc taken from the published papers or GCN reports. 

The afterglow lightcurves of the optical-selected and 
the X-ray-selected samples are presented in Figures 1 
and 2, respectively. The early X-ray afterglow lightcurves 
for the optical-selected sample, and the optical afterglow 
lightcurves for the X-ray selected sample are also present, 
if they are available. In the X-ray selected sample, only 
few cases have simultaneous optical observation. Nine- 
teen out of 20 GRBs in the optical-selected sample have 
simultaneously X-ray observations. We find only one 
case, GRB 080319C, that shows a clear achromatic onset 
bump in both the optical and X-ray bands. This suggests 
that the external shock emission indeed contributes to 
both the optical and the X-ray band. For most optical- 
selected sample GRBs, on the other hand, the early X- 
ray lightcurves either show erratic X-ray flares or inter- 
nal plateaus that are believed to be powered by the GRB 
central engine, or the X-ray observations started only af- 
ter the optical bump peak. Inspecting the two samples 
shown in Figures 1 and 2, we find that the onset bumps 
in the optical-selected sample are usually smoother than 
those observed in the X-ray-selected sample. Seventeen 
out of the 20 GRBs in the optical-selected sample and 
two out of 12 GRBs in the X-ray-selected sample have 
rcdshift measurements. In the following, we mostly use 
only the z-known optical-selected sample in our analysis 
but will use the z-known X-ray-selected sample to con- 
firm the findings. 

We fit the lightcurves with an empirical model pro- 
posed by Kocevski & Liang (2001), 

where F p is the maximum flux at t p , io is a reference 
time, and r and d are the rising and decaying power- 
law indices, respectively. An IDL routine called mpfit- 
fun.pro is employed for our fitting. This routine performs 
Levenberg-Marquardt least-square fit to the data for a 
given model. It optimizes the model parameters so that 
the sum of the squares of the deviations between the data 
and the model becomes minimal. The time interval and 
the fitting curve for each GRB are shown in Figures 1 and 
2, and the fitting parameters are summarized in Tables 1 
and 2 4 . Notice that the lightcurves of some GRBs, such 
as 050820A, 060607A, 070411, 071031, 080330, show sig- 
nificant energy injection or re-brightening features after 
the deceleration bump. We make our fits only around 
the bump. Our fits are also shown in Figures 1 and 2. 

3. CHARACTERISTICS OF THE ONSET BUMP AND 
THEIR CORRELATIONS 

For each fitting lightcurve, we take the full-width-at- 
half-maximum (FWHM) as a characteristic width («;) 
of the bump, and measure the rising and decaying 

4 The reduced x 2 of our fits for some optical lightcurves are 
large. This is due to the fluctuations in the lightcurves and small 
observational errors in the optical data. 
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timescales (t r and t d ) at FWHM. We also derive the ra- 
tios of t r /tp and tr/td- The results are reported in Tables 
1 and 2. The distributions of r, d, t p ,t r , t d , F p , w, t r /t p , 
and t r /td are shown in Figure 3. 

The following statistics applies to the optical-selected 
sample. As shown in Figure 3, the rising index r of most 
bursts are in the range of 1 — 2, with three exceptional 
cases, i.e., GRBs 080330 (r = 0.34 ± 0.03), 060607A 
(r = 4.15 ±0.22), and 050820A (4.45 ±0.76). Inspecting 
the optical lightcurve of GRB 080330, the optical flux 
increases slowly, keeping almost a constant in 300 — 1000 
seconds post the GRB trigger. This feature is similar to 
that observed in GRB 060614. For GRBs 060607A, and 
050820A, their optical lightcurves show a similar behav- 
ior, with rapid increase prior to the peak and a normal 
decay post the peak. The decaying index d is distributed 
in the range of 0.44— 1.77, with a mean value 1.16±0.34. 
Except for GRBs 080330 and 061007, the decay indices 
are well consistent with the isotropic forward shock mod- 
els in a constant density medium. The decay indices of 
the two exceptions are ~ 1.7, slightly steeper than the 
normal decay slope predicted by the isotropic forward 
shock models. The t p is in the range of 10 2 — 10 3 seconds 
with a median value ~ 380 seconds. The distribution of 



F p ranges in 10 



-13 



10 8 ergs cm 2 s 1 , with a mean 



7.25 x 10~ 12 erg cm -2 s _1 . The width w is distributed in 
10 2 — 10 3 seconds. The t r and td peak around 10 2 seconds 
and 10 3 seconds, respectively. The ratio t r /td is narrowly 
distributed in 0.1 — 0.3. However, the distribution of the 
ratio U/tp is much wider, ranging in 0.3 — 1. 

We show various correlations among the deceleration 
parameters of the optical selected sample in Fig. 4, and 
summarize the linear correlation coefficients from pair 
Spearman correlation analysis in Table 3. Tight correla- 
tions are found among t r , td, t p and w. The coefficients 
of these correlations are larger than 0.93. These correla- 
tions are 

\ogt d = (0.48 ± 0.13) + (1.06 ± 0.06) logt r (2) 
\ogt d = (-0.09 ± 0.29) + (1.17 ± 0.11) \ogt p (3) 
\ogt r = (-0.54 ± 0.22) + (1.11 ± 0.08) logt p (4) 
\ogw = (0.05 ± 0.27) + (1.16 ± 0.10) logt p (5) 
\ogw = (0.61 ± 0.11) + (1.05 ± 0.05) logt r (6) 
logw = (0.15 ± 0.02) + (0.98 ± 0.01) logi d . (7) 
The tight t r — td correlations suggest that the struc- 
ture of the bumps among bursts are similar, indicating 
a universal physical origin. No correlation of the decay 
index d with the other parameters is found. This is con- 
sistent with the expectation of blastwave model, in which 
the decay slope is dictated by the density profile and the 
electron spectral index p but has nothing to do with the 
afterglow onset details. On the other hand, the rising in- 
dex r is tightly anti-correlated with both the ratio t r /t d 
and tr/tp, although it is not correlated with t r and t d - 
These correlations read 

log r = (-0.21 ± 0.06) - (1.68 ± 0.19) \ogt r /t p (8) 
logr = (-0.15 ± 0.02) - (0.48 ± 0.05) logt r /t d . (9) 
The w — t p correlation indicates that the wider bumps 
tend to peak at a later time. In addition, both w and 
tp in the burst frame are anti-correlated with the peak 
luminosity at the bump peak L p , 

logL P)0 = (54.6 ± 0.8) - (2.48 ± 0.38) \ogt' p (10) 



log£ p ,o = (0.82 ± 0.79) - (2.16 ± 0.31) logw . (11) 

These results suggest that a dimmer bump tends to peak 
at a later time with a longer duration. 

4. INITIAL LORENTZ FACTOR CONSTRAINTS AND THE 
To - £ 7 ,iso CORRECTION 

The observation properties of the early optical bumps 
seem to be consistent with being due to the onset of the 
external shock afterglow in the thin shell regime. To test 
this hypothesis, we check the correlations of L Pt o and 
the cosmic-frame peak time t p _ z — t p /(l + z) with the 
isotropic gamma-ray energy E 7 ,[ so in Figure 5. We find 
that they are strongly correlated, i.e. 

lo gj L p ,o = (47.10 ± 0.12) + (1.17 ± 0.13) log£ 7 , iso Xi2) 
\ogt p , z = (2.25 ± 0.07) - (0.38 ± 0.07) log£ 7 , iso Xi3) 

with a correlation coefficient k = 0.88 and chance prob- 
ability < 10~ 4 for both correlations. We also derive 
the isotropic optical energy release from 10 seconds to 
10 5 seconds post the GRB trigger (.Eo,iso) using the fit- 
ting curves, and show the correlation between i^o.iso and 
-E 7 ,iso in Figure 5. This correlation, i.e. 

log E v , = (49.63 ± 0.09) + (0.74 ± 0.10) log£ isc , 52 (14) 

has a correlation coefficient k — 0.89 and chance proba- 
bility < 10~ 4 . 

These tight correlations indicate that a GRB with a 
larger E lt i so tends to have a brighter optical afterglow 
peaking at an earlier time. This is well consistent with 
the afterglow onset theory. The shape of the lightcurve 
is consistent with the thin shell case (cf. the thick shell 
case, see Kobayashi et al. 1999; Kobayashi & Zhang 
2007). We therefore apply the standard afterglow model 
in a constant density medium (Sari & Piran 1999) to 
derive the initial Lorentz factor 



r n = 2 



3E~, 



1/8 



32Tmm p c 5 r)t3 
and the deceleration radius 



~ 193(n?7) 



-1/8. 



E. 



7,iso,52 | 

(15) 



R dcc = 2ct p T z p /(l + z) = 2.25 x 10 lb cm r 2 2 t p , z , 2 , (16) 

where rj — E 1 ,- mo / Ek,\so is the ratio between the isotropic 
gamma-ray energy and the isotropic blastwave kinetic 
energy. The results are rather insensitive to n and rj. In 
the following analysis, we take n — 1 cm~ 3 and r\ = 0.2 
throughout. 

With the data reported in Table 1 , we calculate r and 
i?dcc for the GRBs in the optical-selected sample. The 
results also reported in Table 4. The distributions of the 
derived To and i?dcc are displayed in Fig. 6. We can 
see that r is typically a few hundreds, and the typical 
deceleration radius is i?dcc = 1 x 10 17 . From Eq. 15, we 
find To depends on both t PtZ and Ei so . As t PiZ is tightly 
correlated with E- ISOi one expects tight relations between 
r and i? 7 ,iso and between r and t p , z . We show the two 
relations in Figure 7. The best fits give 

log T = (3.69 ± 0.09) - (0.63 ± 0.04) log i P:Z (17) 

with the correlation coefficient n = —0.97 (chance prob- 
ability p < 10~ 4 ) and 

logr = (2.291 ± 0.002) + (0.269 ± 0.002) log£; 7 iso 52 

(18) 
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with correlation coefficient n — 0.91 (chance probabil- 
ity < 1CP 4 ). These tight correlations suggest that t V)Z 
and i? 7! i S o are good indicators of T (hence i?dec)- In 
particular, the latter correlation can be translated to 

r ~ i95£££, iB2 , (19) 

which can be very useful to understand GRB physics (see 
below) . 

5. DISCUSSION 

5.1. Implications for T measurements and the 
T — E 7 ^ so correlation 

Using a sample of GRBs that show the afterglow onset 
feature in the early optical/X-ray afterglow lightcurves, 
we manage to constrain T for a good sample of GRBs 
which can be used to perform a statistical study of Tq for 
the first time. Using a different method (the opacity con- 
straint), the Fermi team recently sets the lower limits of 
r for a number of bright LAT GRBs, e.g. T min ~ 800 
for GRB 080916C (Abdo et al. 2009a), T min - 1200 for 
GRB 090510 (Abdo et al. 2009b), and T min ~ 1000 for 
GRB 090912B (Abdo et al. 2009c). For GRB 080916C, 
one has £ 7iiso ~ 8.8 x 10 54 erg (Abdo et al. 2009a), which 
corresponds to r ~ 1210 according to the r — -E 7j ; so re- 
lation (Eq.[19]). This is consistent with the lower limit 
derived by Abdo et al. (2009a). For GRB 090902B, one 
has £ 7jiso ~ 3.63 x 10 54 erg (Abdo et al. 2009c), which 
corresponds to r ~ 960, which is smaller than r 0! mi n 
derived by the Fermi team. However, as discussed in 
the introduction, we believe that this is due to that the 
Fermi team did not use the most conservative method to 
derive To^m- Since GRB 090902B clearly shows a dis- 
tinct non-thermal component extending to high energy 
which clearly has a different origin from the MeV com- 
ponent (and very likely from different emission regions, 
A. Pe'er et al. 2010, in preparation). One should have 
used the maximum photon energy of the MeV compo- 
nent to estimate To^mm, which would lead to a consistent 
result with the prediction of the r — i? 7 j S o relation. For 
the LAT short GRB 090510, one has E yAso = 3.5 x 10 52 , 
which corresponds to Tq ~ 270 from the To — E^^ so rela- 
tion. This is significantly smaller than To^mm = 1200 set 
by the 31 GeV photon detected during the first second. 
Assuming an external shock origin of the GeV emission, 
Ghirlanda et al. (2009) even derived T ~ 2000 for this 
burst. The inconsistency is significant. This can be again 
due to the non-conservative approach of the Fermi team 
(since the 31 GeV photon is from a different component), 
but more probably it could be that short GRBs do not 
satisfy the r — E 1 - lso relation derived for long GRBs. 

With T = 100 - 1200 derived from the r - E lsor/ 
relation, we expect that the corresponding observed t p is 
in the range of 30 ~ 1400 seconds for a typical redshift 
z = 2 according to the r — t p , z correlation (Eq.[17]). 
In this time period, the observed X-rays are generally 
dominated by the GRB tail emission or flares. This may 
be, at least partially, the reason why not many early 
X-ray lightcurves show the clear afterglow onset bump 
signature. 

The T — -E 7 ,iso relation is very useful to pin down the 
prompt emission physics of GRBs. One interesting ob- 
servational correlation is the Amati relation, i.e. E p oc 
E* iso (or E p oc L* iso , with k ~ (0.4 - 0.5), both as a 



bulk correlation among bursts and an internal correlation 
within a burst (Amati et al. 2002; Wei et al. 2003; Liang 
et al. 2004; Yonetoku et al. 2004; Lu & Liang 2009). 
However, all the prompt GRB emission models predict 
E p as a function of both i? 7l i S o (or L^s so ) and r (e.g. 
Table 1 of Zhang & Meszaros 2002). There is no straight- 
forward theory that predicts the relationship between Tq 
and i? 7 .iso (or L^j so ). As a result, any theoretical model 
can be argued to interpret the Amati relation, given a 
designed r — E 7 ^ so correlation. The r — -E 7l i so corre- 
lation discovered here therefore poses great constraints 
on many prompt emission models. For example, the in- 
ternal shock synchrotron model predicts E p oc L 1 / 2 Fq 2 
(Zhang & Meszaros 2002). The Amati relation essen- 
tially requires that To oc L°. The E p oc E®-?J a relation, 
combined with the trivial proportionality L oc E 7 ^ so (a 
non-correlation between GRB duration and luminosity) , 

would lead to E p oc E^E^ 4 oc £~° s ° 4 , which means 
that Ep is essentially constant for different -B 7 .iso values. 
This is in contradiction with the Amati relation, which 
can be regarded as another argument against the inter- 
nal shock synchrotron emission model of GRB prompt 
emission (see also Kumar & McMahon 2008; Zhang & 
Pe'er 2009). 

5.2. Early Optical vs. X-ray emission: different physical 

origins? 

The simultaneous observations in the optical and X- 
ray bands during the early afterglow phase also hold the 
key to address whether the broad band emission is from 
the same emission component. With the prompt slew- 
ing capability, the X-ray Telescope (XRT, Burrows et 
al. 2004) on board Swift has established a large sam- 
ple of X-ray afterglow of GRBs. Generally, the XRT 
lightcurves are composed of a few power-law decaying 
segments and some erratic flares. Although the X-ray 
lightcurves are diverse among bursts, they can be roughly 
classified into three groups with a large, uniform sample 
established by XRT. The majority is the so-called canon- 
ical XRT lightcurves characterized by a steep-shallow- 
normal-steep decay pattern (Zhang et al. 2006; Nousek 
et al. 2006; O'Brien et al. 2006), although not all the 
segments show up in every burst (Evans et al. 2009). The 
second group is composed of those lightcurves that show 
a single power-law decay from early to late epochs (Liang 
et al. 2009; Evans et al. 2009). The third group includes 
some GRBs that show an "internal plateau" that is fol- 
lowed by a rapid drop with a decay slope steeped than -3 
(Liang et al. 2007; Troja et al. 2007; Lyons et al. 2009). 
The physics origin of the X-ray emission is still a mystery 
(e.g. Zhang 2007 for a review), although the consensus 
is that there might be diverse origins (Liang et al. 2007) . 
It is clear that the X-ray flares and internal plateaus are 
of an internal origin. However the origin of the canonical 
lightcurve is still subject to debate. "Closure" -relation 
analyses suggest that the normal decay segment following 
the shallow decay one in the canonical XRT lightcurves 
are roughly consistent with the forward shock models 
(Willingale et al. 2007; Liang et al. 2007), favoring the 
long lasting energy injection models for the shallow-decay 
segment. However, the optical/X-ray chromatic behav- 
ior around the shallow and normal decay transition time 
(Panaitescu et al. 2006; Fan & Piran 2006; Liang et al. 
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2007) suggests that the X-ray and optical emissions are 
two independent components. Some models attribute 
the entire X-ray emission to the late emission from the 
central engine, probably related to the long-term accre- 
tion history of the central engine (Ghisellini et al. 2007; 
Kumar et al. 2008; Cui et al. 2009; Cannizzo & Gehrels 
2009), but the consistency with the "closure-relation" of 
the external shock model predictions is not naturally ex- 
plained in these models. Interestingly, Yamazaki (2009) 
recently suggested that the X-ray emission might be an 
independent component prior to the GRB trigger, and 
that the apparent shallow-to-normal transition is merely 
reference time effect. Liang et al. (2009) systematically 
studied the Swift canonical GRB lightcurves and con- 
firm that shifting the reference time can indeed stretch 
the canonical lightcurves to single power law lightcurves, 
and proposed a unified picture for the physical origin 
of both the canonical and the single power-law decaying 
XRT lightcurves. 

As shown in Figure 1, the optical and X-ray lightcurve 
behaviors are dramatically different at the early epoch, 
i.e., t < 1000 seconds post the GRB trigger. This 
strongly suggests that the radiations from the two energy 
bands are not from the same component. The tight cor- 
relations between the observables of the prompt gamma- 
rays and early optical afterglows shown in Fig. 5 strongly 
suggest that the optical emission is likely the "afterglow" 
of the GRB fireball (external shock component). This is 
consistent with the smooth afterglow onset feature ob- 
served in the optical band for these GRBs. On the other 
hand, it also suggests that the early bright X-ray emis- 
sion is not from the external shock. One natural question 
would be: where is the external shock X-ray component? 
Inspecting the details of the X-ray and optical afterglow 
lightcurves in Figure 1, we can see that this component is 
very likely hidden underneath some brighter X-ray emis- 
sion components in the early epochs (e.g. flares, internal 
plateaus, or even normal plateaus). Interestingly, one 
can find an X-ray decay slope similar to that of optical in 
the late epochs in half of GRBs in our sample, including 
GRBs 050820A, 060418, 060605, 061007, 070318, 070411, 
071031, 080319C, 080810, 081203, and 090102. We there- 
fore cannot exclude the possibility that the late X-ray 
and optical emissions share the same external shock ori- 
gin. 

The mixture of different emission components in the 
X-ray observations (Willingale et al. 2007; Liang et al. 
2007; Liang et al. 2009; Nardini et al. 2009) also natu- 
rally interprets the fact that the rarity of the onset after- 
glow feature observed in the X-ray band. This requires 
that other early X-ray components are not bright enough 
to outshine the external shock component. In fact, we 
find only 13 out of ~ 400 cases in the current XRT 
lightcurve sample. They are shown in Figure 2. Since 
the optical band is less affected by the other emission 
components related to the central engine, one naively 
expects that the X-ray onset cases should have achro- 
matic optical onset feature as well. Unfortunately, the 
optical/X-ray joint observations in this sample are rare: 
only GRB 080507 and GRB 080319C have early optical 
observations. The optical lightcurve of GRB 080507 is 
sparse. For GRB 080319C one indeed observes an achro- 
matic onset feature. There is an earlier decay feature in 
the optical lightcurve of GRB 080319C. It may be asso- 
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ciatcd with an internal emission component within such 
an interpretation. 

6. CONCLUSIONS 

We have made an extensively search for the afterglow 
onset "bump" feature from early afterglow lightcurves, 
both in the optical band (through literature survey) and 
in the X-ray band (through systematically analyzing the 
Swift XRT data). Twenty GRBs are identified in the 
optical-selected sample and 12 GRBs are found in the 
X-ray-selected sample. We fit the onset bumps with a 
smooth broken power-law and measure their characteris- 
tics. The rising index r for most bursts is 1 — 2, and the 
decay index d is 0.44 — 1.77. These are well consistent 
with the forward shock models. The peak time t p is in 
10 2 — 10 3 seconds with a median value of ~ 380 seconds. 
The width of the bumps measured at FWHM is 10 2 - 10 3 
seconds, and the typical rising time t r and decaying time 
td are 10 2 seconds and 10 3 seconds, respectively. The ra- 
tio of t r /td is narrowly distribution around 0.1 — 0.3, and 
the ratio t r /t p has a distribution in the range of 0.3 — 1. 

Most GRBs in our optical-selected sample have red- 
shift measurements. We analyze pair correlations among 
the bump characteristics. We find that the pulse width, 
rising time, decaying time, and the peak time are 
strongly correlated. Bumps that peak later are dimmer 
and wider. No correlation between the decay index d 
with other parameters is found, but the rising index r 
is tightly anti-correlated with both the ratio t r /td and 
t r /tp, although it is not correlated with t r and td- 

We analyze the relation of the optical afterglow bumps 
with prompt gamma-ray properties. We find that a GRB 
with larger -E 7i i so tends to have a brighter optical after- 
glow, and tends to be decelerated by the surrounding 
medium earlier. These tight correlations strongly sug- 
gests an external shock afterglow origin of the early op- 
tical emission. 

Within the framework of the standard forward model 
in a constant density circumburst medium, we calculate 
the initial Lorentz factor r and the deceleration radius 
i?doc for the GRBs in the optical-selected sample. The 
derived Tq ranges from 100 to about 600, while i?dcc 
narrowly distributed around 10 cm. 

Intriguingly, we discover a tight correlation between Tq 
and -E 7 ,iso- For typical values n = 1 cm~ 3 and rj = 0.2 
(r is insensitive to the values of n and rj, we obtain the 
correlation Eq.(19). This correlation is very important 
to understand GRB prompt emission physics, and may 
serve as an indicator of Tq for other long GRBs. In par- 
ticular, the correlation disfavors the internal shock syn- 
chrotron emission model for the GRB prompt emission. 
The extrapolation of the correlation is consistent with 
the r , min of GRB 080916C derived by the Fermi team 
(Abdo et al. 2009a), but are inconsistent with the To m j n 
of GRB 090510 and GRB 090902B derived by the Fermi 
team (Abdo et al. 2009b, c). We point out that this is 
because the Fermi team did not use a more conservative 
approach to set ro jm i n . 

The X-ray-sclcctcd sample only has two cases with red- 
shift measurements. The derived r 's from these two 
cases are also consistent with the r — i? 7 ,i so correlation 
derived from the optical-selected sample. There is one 
case (GRB 080319C) that shows an achromatic afterglow 
onset feature. 
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Most optical-selected sample has early X-ray emission 
components not from the external shock. This reinforces 
the diverse origin of early X-ray afterglow in most GRBs. 
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XRT observations and fitting result for our sample 
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TABLE 3 

The optical and X-ray with know redshift in our sample 
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TABLE 4 

Spearman correlation coefficients of the Optical Selected Sample 



L Pi o t p r a t r t d t r /t d tr/tp w 



L Pi0 -0.90 X X -0.89 -0.88 X X -0.88 

tp X X 0.95 0.93 X X 0.94 

r XXX 0.90 0.90 X 

d X X X X X 

t r ' 0.98 X X 0.98 

X X ~1 

tr/td X X 
tr/tp X 



GRB initial Lorcntz factor and Tq — -E 7 .i SO relation 



9 




Fig. 1. — Lightcurves of the optical-selected sample (dots). Lines are the best fits with a smooth broken power-law to the optical data 
in a selected interval marked as a red solid segment. Blue dashed segments just extend the fitting curves in a broad time regime. The 
simultaneous X-ray data observed with Swift/XRT (crosses with error bands) are also present. 
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Fig. 1 — continued 
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FlG. 2. — The same as Figure 1, but for the X-ray selected sample. 
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Fig. 3. — Distributions of the characteristics of the optical-selected sample. 




Fig. 4. — Pair correlations among the characteristics of the optical-selected sample. Lines are the best fits. 
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FlG. 6. — Distributions of Tq and R d for the optical-selected sample 



GRB initial Lorentz factor and Tq — -E 7 ,iso relation 
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Fig. 7. — Relations of To to £^ so T , ^p,0> an d + 2). Lines are the best fits for the optical-selected sample. GRBs 070208 and 

0803f9C in the X-ray selected sample have redshift measurement. We mark the two GRBs with red stars. 



